Detecting light long-lived particle produced by cosmic ray 
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We investigate the possibility of detecting light long-lived particle (LLP) produced by high energy 
cosmic ray colliding with atmosphere. The LLP may penetrate the atmosphere and decay into a 
pair of muons near/in the neutrino telescope. Such muons can be treated as the detectable signal for 
neutrino telescope. This study is motivated by recent cosmic electron/positron observations which 
suggest the existence of O(TeV) dark matter and new light O(GeV) particle. It indicates that dark 
sector may be complicated, and there may exist more than one light particles, for example the dark 
gauge boson A' and associated dark Higgs boson h' . In this work, we discuss the scenario with 
A' heavier than h! and h! is treated as LLP. Based on our numerical estimation, we find that the 
large volume neutrino telescope IceCube has the capacity to observe several tens of di-muon events 
per year for favorable parameters if the decay length of LLP can be comparable with the depth of 
atmosphere. The challenge here is how to suppress the muon backgrounds induced by cosmic rays 
and atmospheric neutrinos. 

PACS numbers: 12.60.-i, 13.85.Tp 



I. INTRODUCTION 

The pursuit of new light boson (LB) with mass much 
lighter than W and Z bosons in the standard model(SM) 
has a long history. A well known example is a light 
gauge boson under an extra U(l) gauge group beyond 
the SM gauge group SU(3) C ® SU(2) L ®U(1) Y . However 
such kind of new LB is stringently constrained since the 
current precise measurements are in excellent agreement 
with the SM predictions. Obviously the interaction be- 
tween the new LB and SM sector should be tiny to evade 
current constraints. With the same reason detecting LB 
experimentally is a very challenging task. 

Recently the LB attracts more attention due to the 
new cosmic observations by PAMELA [J, ATIC and 
Fermi Q . PAMELA collaboration reported excess in the 
positron fraction from 10 to about 100 GeV but absence 
of anti-proton excess pj. It is still consistent with the 
results released by ATIC @ or Fermi Q experiments. 
The afterwards investigations suggested that these new 
observations need new source of electron/positron which 
may come from the dark matter (DM). The DM parti- 
cles might annihilate or decay into electron/positron in 
the halo today. Moreover man began to realize the pos- 
sible connection between the heavier DM at O(TcV) and 
the LB at O(GeV). If the annihilating DM produces the 
extra electron/positron, there is a mismatch between the 
DM annihilation cross section expected in the epoch of 
freeze-out and that required to account for recent new 
observations. Namely the present DM annihilation cross 
section is too small. The new O(GeV) LB, via the so- 
called Sommerfeld enhancement, can fill the gap. More- 
over, in order to be consistent with only electron/positron 
excess, LB decays preferably into charged leptons [||. 

Generally speaking such new boson may be scalar, 
pscudoscalar or gauge boson. The interactions between 



the LB and SM sector could arise from the mixing be- 
tween the LB and photon or Higgs. It is quite interest- 
ing to investigate how to detect such kind of LB. Many 
authors have studied how to produce such LB and de- 
tect them at colliders, namely low-energy collider, large 
hadron collider or fix-target experiment. For the collider 
search, one often requires the life time of LB is short, as 
a consequence the charged leptons as the LB decay prod- 
ucts could be observed at the detector. Such charged lep- 
tons could be clearly and easily identified. The construc- 
tion of realistic model of dark sector showed that dark 
sector can be more complicated. The assumption that 
LB is short-lived might be incorrect. If dark sector con- 
tains an array of LBs besides a light gauge boson, some 
lighter ones may be long-lived due to the suppressed in- 
teraction with SM sector. Several recent works provided 
an interesting approach to search such long-lived parti- 
cle (LLP) [111. The DM trapped inside the Sun/Earth 
would annihilate into LBs. If the LB has a long life- 
time, it can travel through the Sun/Earth and decay into 
gamma rays or charged leptons which could be observed. 

In this paper we point out another possible way to 
search such kind of LLP via the high energy cosmic rays. 
The high energy cosmic rays interact with atmospheric 
nucleons every second and can be treated as a natural 
and costless high energy hadron collider. If a proton 
with energy of E ~ 10 4 GeV in cosmic ray collides with 
an atmosphere nucleon, the center-of-mass energy y/s is 
approximately ^/2m^E ~ 10 2 GeV. Such mechanism 
can copiously produce LLPs through pN collisions. If 
the lifetime of LLP is appropriate, it can penetrate the 
atmosphere and arrive at the neutrino telescope. The 
subsequently decay of LLP into a pair of muons can be 
observed by the telescope [1, |T(| ■ 

If LLP decays near the detector, the difficult task is 
to distinguish the signal muon pair from the huge muon 
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backgrounds which are produced by high-energy cosmic 
rays through hadron(for example 7f) decays and/or QED 
process. Provided that the time information of muon 
is precisely recorded, one can identify two signal muons 
which arc supposed to arrive at detector at the same time, 
not heavily polluted by two irrelevant coincident paral- 
lel muon events. In order to suppress the muon back- 
grounds, the analysis focused on quasi-horizontal events 
might be important [U, [HI . A more optimistic case is 
that LLP decays inside the detector with an 'obvious' 
decay vertex. The challenge here is how to distinguish 
almost parallel di-muon from the single muon. Another 
possible case is that for a single collision between cosmic 
ray particle and atmosphere nucleon, more than one LLP 
are produced. In this case the multi muon pairs appear 
in detector at the same time will be a significant signal. 

In this work, we assume the LB has a long lifetime 
to penetrate through the atmosphere, some mechanisms 
which satisfy this requirement will be discussed. We con- 
sider the LB production by the cosmic rays and simulate 
di-muon events from LB decay. In order to detect such 
high energy muons, we focus on the large volume neu- 
trino telescope IceCube [l3| which will reach an effective 
detecting area in square kilometers. The detector are in- 
stalled under the ice surface in a depth of 1.4 km in order 
to suppress muon backgrounds produced in atmosphere. 
Moreover, there is a extension of IceCube namely Deep- 
Core which is still under construction (l4j . DeepCorc will 
be installed in the more deeper location which will suffer 
less above-mentioned muon backgrounds. 

This paper is organized as following. In the section II, 
we describe two classes of models which contain LLPs, 
discuss the main LLP production processes and calculate 
the LLP production cross section in pN collision. We 
utilize PYTHIA to do the calculation and simulate the 
LLP events. In Section II, we also investigate the LLP 
production flux produced by primary cosmic ray. We find 
that for a LLP production process with O(10) pb, the 
production rate of LLP can reach to O(10 4 ) per square 
kilometer and per year. In Section III, we investigate the 
possibility to detect di-muon signal from LLP decay in 
the neutrino detector. The conclusions and discussions 
are given in the last section. 



II. THE PRODUCTION OF LLP 
A. The model 

In the models discussed in Ref.[^,Q, the dark sector in- 
cludes both weakly-interacting-massive-particle (WIMP) 
and LB under a certain new gauge group. For the sim- 
plest case with an extra U(l) group, the dark sector can 
interact with SM sector though kinetic mixing, namely 
a new U(l) gauge field A'^ could mix with U(l)y field 
in the SM. In addition there is another possible mix- 
ing between the SM Higgs field H and scalar h! which 
will break extra U(l) group to induce mass to A' . The 



Lagrangian can be written as, 

c = -^ 2 + ^ + MI 2 - V(h') 

+ \ h , H {h^h'){H^H)+C DM + CsM. (1) 

Here k and \h'H are the two mixing parameters which 
will be determined (constrained) by experiments. Cum 
is the lagrangian of DM which includes the 0(TeV) DM 
kinetic and mass terms and its gauge interaction. 

After the spontaneous breaking of extra U(l) group, 
the A' and hi will get the mass of tua' and ■ If 
scalar mixing Xh'H is neglected, the lifetime of LB 
is determined by and gauge kinetic mixing pa- 

rameter k. For light gauge boson mainly decay into 
charged leptons, the travel distance can be approxi- 
mately estimated as I = jct ~ 7c/ ' {an 2 mA') ~ 10~ 5 m 
(7/10 3 )(K/10- 3 )- 2 (m A VlGey)- 1 . Typically the light 
gauge boson is not a long-lived particle. However for 
the dark Higgs boson and provided that m' A > m! h , 
dark higgs will decay into SM fermions through the tri- 
angle diagrams. The decay width is suppressed by a 
factor of k~ 4 @. Therefore the dark higgs can be a 
typical long-lived particle (LLP) with the decay length 
as large as 10 7 km ||. In such kind of models, the 
main LLP production mode is scalar-strahlung process 
pp — s- A'* + A — > A'h' + X ( X denotes anything). 

The dark sector could be more complicated. As dis- 
cussed in Ref. [HI, if the dark sector has a more 
complex gauge group configuration, there might exit a 
series of LBs including the LLP. In such case, the high 
energy pN collisions may copiously produce the extra 
gauge boson A' and the production rate depends on the 
mixing parameter k and tua'- A' subsequently decays 
quickly into LLP, namely the dark Higgs boson h! . The 
process can be depicted as pp — > A' + X — > hi o! + X 
(a' represents another light gauge boson or light pseudo 
scalar which we do not discuss further its feature here. 
For simplicity we take the mass of hi and a! to be equal). 
If LLP propagates some distance which depends on its 
lifetime and decays into charged leptons as motivated by 
cosmic electron/positron data, we can utilize neutrino 
telescope to observe such kind of muons. Note that the 
LLP interacts with usual matter weakly, it can penetrate 
into the Earth without loss of energy 

In this case, we do the calculation as model- 
independent as possible and the input parameters are 
chosen as kinetic mixing k and the mass of extra gauge 
boson m,A' which determine A production rate, mass 
nth' and lifetime r of h'. Typically branching ratio of A' 
into dark sector is much bigger than that into SM sector. 
Thus in our numerical simulation we assume that A' de- 
cays only into hi . We also assume the branching ratio of 
hi into muon pair is 1. In fact, such branching ratios are 
calculable in the model we mentioned above, the detailed 
calculations can be found in the Ref. @, Q . Here we just 
treat them as free parameters, and our results could be 
adjusted to satisfy specific model by multiplying a factor 
of Br(A' -> h')Br(h' -> 
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B. Simulation of LLP production 



and set the factorization scale Q 2 = s. To do a Monte 
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FIG. 1: Cross section of LLP production in pp collision as a 
function of cosmic primary proton energy. The dash lines and 
the dashed-dot lines represent process qq — > A' and qq —¥ h'A' 
respectively. For the lines with same shape, the upper and 
lower one represent three benchmark point (m A i, m h i) — (1.2 
GeV, 0.4 GeV), (2.1 GeV, 0.7 GeV), (3.0 GeV, 1.0 GeV) 
respectively. 

We utilize PYTHIA [l7| to simulate the LLP produc- 
tion events and calculate the cross section. The main 
free parameters are the mass of A' and hi , kinetic mixing 
parameter k, and dark sector gauge coupling a' for the 
Higgs-stahlung process. Here we choose three benchmark 
points with {m A >,m h >) = (1.2 GeV, 0.4 GeV), (2.1 GeV, 
0.7 GeV), (3.0 GeV, 1.0 GeV) and a' = a for simplic- 
ity. The k is severely constrained by experiments, and we 
choose k as 10~ 3 which is still allowed [Toj . In addition, 
for the case < itla> we discussed here, the k could 
not be very small to avoid destroying the success of BBN. 
The lifetime of h! must be shorter than 1 second, and the 
k should not be very smaller than 10 -4 @. 

For the scalar-strahlung process, the cross section at 
parton level is given by (neglecting the mass of parton), 



x 2 are randomly 

f/s < Xl < 1, 



Carlo calculation [L9jj , the variables x\, 
chosen within the ranges (m^ + 

(rn,A> + TUh') 2 1 '(xis) < X2 < 1. In order to improve the 
convergence, we technically define new integration vari- 
ables as dV = dhixidhix 2 . The integration is then 
transformed into 
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Xif^{xi)x 2 f q N {x 2 )d{s) q qdV' 



(4) 

where i denotes i-th configuration of cross section. 

For the single A' production which is a 2 — > 1 pro- 
cess, the cross section contains one 5 function which fixes 
s = m 2 A ,. We transform integration variables X\, x 2 to 
s, y through dx\dx 2 = dsdy/s and integrate out the 5 
function. The total cross section can be written as 
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xifp{xi)x 2 f q N {x 2 )dy (5) 



where Xi = m A ie v j\fs and x 2 = rriA'e v j\fs. From 
^1,2 < 1, we choose y in the range of — \n(yfs / m a>) < 
y < ln(y/s/m A >). 

We show the LLP production cross sections for the pp 
collision in Fig. [TJ From the figures we can see that the 
scalar-strahlung process cross sections are much smaller 
than those of A' resonance process. It is simply because 
the scalar-strahlung process is the 2 — > 2 process which 
has more power of coupling and smaller phase space. Nu- 
merically in the energy region of O(100) GeV, the single 
A' production is several orders of magnitude larger than 
that the scalar-strahlung process. From the Fig. [1] we 
can also conclude that cross section is very sensitive to 
the mass of A'. This is quite understandable provided 
the quick rise of qq luminosity for the lower mass A'. 
Note that flux of cosmic ray decreases quickly with the 
increment of its energy, thus the behavior for low energy 
cosmic ray plays the major role to produce LLP. 

In our simulations, we randomly generate muon pairs 
in the h! rest frame, then boost them to the A! rest frame 
and lab frame. 
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where s = 
tons with 



m. 



x\X 2 s is the center of mass energy of par- 
momentum fractions of x\ and x 2 . k = 
Am 2 A ,m 2 L , I '(2-\/i) is the momentom 
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of A' in the center of mass frame. The total cross section 
is 



(3) 



where dV represents dx\dx 2 , q denotes sum over all the 
quark and anti-quark and / is the parton distribution 
function (PDF). We use the CTEQ6M PDF [3 here 



C. Flux of long-lived particle produced by high 
energy cosmic rays 

In this subsection, we evaluate the LLP production 
rate by high energy cosmic ray. From several GeV to 
energy range above 10 6 GeV, the flux of primary nucleons 
in the cosmic rays is approximately written as [20| 



<f> N (E) f« l.8(E/GeV)- 
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i 2 s sr GeV ' 
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where the differential spectral index a is 2.7 under 10 6 
GeV, 3.0 from 10 6 GeV to 10 10 GeV, and finally 2.7 above 
10 10 GeV. The main component of primary cosmic nu- 
cleons is proton. To produce energetic LLP, the energy 
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FIG. 2: The LLP differential production rate as a function of 
the energy of primary cosmic ray. The notations are same as 

Fig.m 



of cosmic rays is required to be as high as or above 10 2 
GeV. For secondary hadrons induced by such energetic 
cosmic ray, the interaction length may be smaller than 
the decay length in the atmosphere. Thus the collisions 
between the atmosphere nucleons and secondary mesons 
induced by primary cosmic hadrons are important 1 . The 
flux of new LLP particles can be estimated by [ll[ 



E 



dE $ h {E)Vk(E), 



(7) 



where h denotes the primary cosmic hadrons and sec- 
ondary nucleons, pion and kion, 4>h is the flux of hadron, 
and Vh (E) is the probability of producing LLP hi in one 



collision. V\,(E) is approximated as Vfr(E) ss Aafcf 4 /a T 
[111 ], where A ~ 14.6 is the average nucleon number of a 
nuclei in air, is the cross section to produce h', and 
ctt is the total cross section of cosmic hadrons in atmo- 
sphere about O(10 2 ) mb which is approximately parame- 
terized as <j T « C* +C' 1 l \n(E/GeV) + C* \n 2 (E/GeV) 2 . 
Requiring final muons from LLP decay with large energy 
and flux, we choose the cosmic ray energy region from 
10 2 GeV to 10 7 GeV. 



For a simple and conservative estimation, this work does not 
include the effects of secondary hadrons. The influence of sec- 
ondary hadrons to final LB flux can be included by multiplying 
an extra O(l) factor. This factor is determined by the flux of 
secondary hadrons llll . 

In fact, the total cross sections of high energy cosmic hadrons 
are not precisely determined. There arc several models which 
are consistent with the exiting observations roughly, but more 
data are needed to give a more concrete result. The detailed 
discussions about different models and parameter selections can 
be found in Ref. |21| . Here we use the result of VENUS model 
l22l given in Ref. |21| as an acceptable approximation from 10 2 
GeV to 10 7 GeV. For present experiment results and comparisons 
between different theoretical predictions, one can see the Ref. 
[231 and references therein. 



We can firstly estimate the LLP production rate in 
order of magnitude as 
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(8) 

The formula indicates that for the single A' production 
cross section of 10 pb, the cosmic ray will produce at 
least 10 4 LLPs per year and per square kilometer. In 
Fig. [2] we show the results of LLP differential rate as 
a function of primary cosmic ray energy. From the fig- 
ure we can see that the event rate can reach O(10 2 ) to 
O(10 4 ) kmT 2 yr~ x sr~ x . However for the scalar-strahlung 
process, the event rate is very low from O(10 _1 ) to O(10) 
due to the small cross section and we do not discuss this 
process below. 



III. DETECTING LONG-LIVED PARTICLE AT 
HIGH ENERGY NEUTRINO DETECTOR 

In this section we will discuss the possibility of detect- 
ing such LLP at neutrino detector. If the LLP penetrates 
through atmosphere and decays into muons near or in- 
side detector, this signal can be observable at neutrino 
telescope. In the last section, we have evaluated the hi 
production rate. For simplicity, we assume that all hi s 
are produced at the upper atmosphere. Moreover, the 
hi flux can be treated as isotropic due to the isotropy of 
cosmic ray. The long-lived h! interacts with SM particle 
weakly, thus it is safe to neglect the energy loss before it 
decays into muon pair. 

The lifetime of hi is model-dependent, and it is treated 
as a free parameter here. The probability of a particle 
decays between two point is given by @: 
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decay 
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-D/l, 



{l-e- d '\ (9) 



where D is the distance between the production point 
and entry point, d is the distance between two point, 
I = jvt is the decay length of particle. Here D is approx- 
imately 10 ~ 20 km, and it is the depth of atmosphere 
plus distance from detector to horizon; d is the neutrino 
telescope's size (if the detector has the capacity to recog- 
nize LLP decay events nearby, d can be larger). If D is 
far larger than d and I, the decay probability is roughly 
e -D/i. jf ^hg decay length I is far larger than D and d, 
the decay probability is d/l. Therefore, it is difficult to 
observe the LLP decay if the decay length of LLP is too 
large or too small. The most promising case is that the 
decay length I is comparable with D. 

As mentioned in the introduction if LLP decays outside 
the neutrino detector, the atmosphere muon background 
is huge. These muons are generated mainly from sec- 
ondary charged pion and kaon decay. Most muons are 
in lower energy regime due to the relatively large cosmic 
ray flux in this energy region and energy lost in matter. 
The formula of atmosphere muon flux is similar to that 
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of cosmic ray which is approximately given by [2C 



O.UE, 



-2.7 



1 



0.054 



cm? s sr GeV \ 



1AE U 



UbGeV 



1 



1.1E,, cos i 



-). 



850GeV 

(10) 

where 9 is zenith angle (here 9 < 70°), and this formula is 
valid when the probability of the muon decay can be ne- 
glected (i.e. Ef, > 100/cos6» GeV). From the Eq. ([T0]) . we 
can see that in the energy region O(10 2 ) GeV the atmo- 
spheric muon flux is only one order of magnitude smaller 
than primary cosmic ray, while the LLP flux is about ten 
orders smaller. Requiring that two muons arrive at detec- 
tor in a tiny time window could effectively reduce back- 
ground which contains two uncorrelatcd muons. Note 
that lots of SM processes can produce muon pair events 
directly (24|. For example, one single shower contains 
many hadrons, and muons from two hadrons decay may 
be treated as a pair of muons. In addition the electro- 
weak Drell-Yan process can also produce muon pair di- 
rectly. In a word, if the LLP does not decay inside the 
detector, it is very challenging to distinguish signal from 
the atmospheric di-muon backgrounds. 

However there is still hope to detect muon pair from 
LLP decay near the detector. The energy of LLP is typ- 
ically less than 1 TeV. In such energy region, most of 
the atmospheric muons are absorbed in the solid/liquid 
matter. For example, for 10, 10 2 and 10 3 GeV muons, 
the muon range is 0.05, 0.41 and 2.45 km.w.e respectively 
[25| . The neutrino detector is often installed in deep un- 
derground with shield of several km.w.e. rock/water/ice. 
Such shield can prohibit the low-energy atmospheric 
muons to arrive at the detector, especially for direction 
with large zenith angle and even quasi-horizontal direc- 
tion. Detecting LLP decay in these directions is more 
promising. Thus if the LLP decays not far from detec- 
tor, the resulting di-muon may be identified. 

To detect the clean LLP signal, we expect that LLP 
happens to decay inside the detector and we can observe 
a pair of "suddenly" appeared muons. In order to recog- 
nize di-muon event, we require that the energy of muons 
must be above the detector's threshold energy. Moreover, 
because the high energy atmospheric neutrino may pro- 
duce single high energy muon when traveling through the 
detector, we also require that the two tracks of di-muon 
should be identified separately 3 . In our simulations, we 
require that the angle between two muon tracks is greater 
than O(10~ 4 ). 

It needs to mention another background arising from 
atmospheric neutrino, which can induce di-muon events 
inside the detector through inelastic vN collision. vN 
collision can produce muon plus charm hadron and 



the charm hadron will subse quen tly induce another 
muon via semi-leptonic decay |26[. Such two muons 
will be identified as a di-muon event. Such di-muon 
event rate may be larger than the LLP decay due 
to large atmospheric neutrino flux (for example, the 
atmospheric muon neutrino flux around 100 GeV is 
about O(10 _4 )m~ 2 ,s _1 sr~ 1 GeV~ 1 ). However, associ- 
ated hadronic shower in l/N collision could be utilized 
to suppress such backgrounds. 
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FIG. 3: Di-muon event rate as a function of di-muon en- 
ergy. The dash lines and the solid lines correspond to differ- 
ent energy threshold of 3.2 GeV for Super-Kamiokande and 
100 GeV for IceCube respectively. For the lines with same 
shape, the upper and lower one represent three benchmark 
point {m A :,m h ,) =(1.2 GeV, 0.4 GeV), (2.1 GeV, 0.7 GeV), 
(3.0 GeV, 1.0 GeV) respectively. 



■ 1 1 1 1 1 1 1 



■ 1 1 1 1 1 11 1 



■ 1 1 1 1 1 




FIG. 4: Di-muon event rate as a function of separation angle 
between two muons. The notations are same as Fig. [3] 



3 The authors of Ref. § have pointed out that it is possible to 
distinguish di-muon signal from single muon even the separation 
of di-muon is not large enough. They provided two handles to 
identify di-muon signal with energy above critical energy utilizing 
the different characteristics of Cherenkov radiation. 



Now we scrutinize the possibility of detecting LLP at 
the large volume Cherenkov detector IceCube. IceCube 
is a lfcm 3 detector which consists of 4800 digital optical 
modules and is installed between depth of 1450 m and 
2450 m in the south pole [ll]. In Fig. [3] and Fig. H 
we show the di-muons rate as the functions of the di- 
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muon energy and the separation angle. Here we assume 
all the LLPs decay in the detector. We take the each 
muon energy threshold of IceCube to be 50 GeV and the 
separated angle of di-muon greater than 10 -4 . For com- 
parison, we also show the rate of di-muon events with 
energy above 3.2 GeV which corresponds to the thresh- 
old at the Super-Kamiokande detector [27| ■ From Fig. [3] 
we find that di-muon event rate with energy larger than 
100 GeV could be O(10 3 ) per square kilometer, per year. 
From Fig. 2] we find that the separated angle of most di- 
muons is about O(10 -2 ). For the Super-Kamiokande de- 
tector with R = 16.9 m, H = 36.2 m and energy thresh- 
old 1.6 GeV, it is suitable for detecting the low energy 
di-muon events with large flux. However it is difficult to 
distinguish the di-muon events from single muons due to 
the limited volume. 
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FIG. 5: Di-muon rate as a function of the lifetime of LLP, 
where {m A >,m h >) = (1.2 GeV, 0.4 GeV), (2.1 GeV, 0.7 GeV), 
(3.0 GeV, 1.0 GeV) from top to bottom. The size of detector 
are taken as 1 km and the zenith angle is less than 9 < 70°. 
We choose each muon energy threshold of detector to be 50 
GeV and the separated angle of di-muon to be greater than 

io- 4 . 



In Fig[3]and Fig. we did not include effects of LLP 
decay. The Fig. [S] shows the di-muon rate as a function 
of the lifetime of LLP. We assume that the LLP is ran- 
domly produced between the altitude of 15 km and 20 
km and the LLP decay length is taken as I = crE/mh'- 
The size of IceCube is chosen as 1 km. The effective area 
of detector is not taken into account since it needs the de- 
tailed simulation including realistic experimental trigger 
and condition. 

The detector is insensitive to the direction for low 
zenith angle. The reason is that for a given lifetime 
t, the LLP has similar I and P decay in these direc- 
tions. Therefore we show the rate in the directions 
with zenith angle below 70°. We choose each muon en- 
ergy threshold of detector to be 50 GeV and the sep- 
arated angle of di-muon to be greater than 10 -4 here, 
and find that the di-muon rate can reach several tens 
for the most optimistic cases. For example, the event 



rate can be 42, 9, 3 per year for three benchmark 
points {m A ',m h >) = (1.2 GeV, 0.4 GeV), (2.1 GeV, 
0.7 GeV), (3.0 GeV, 1.0 GeV) if the lifetime of LLP is 
~ 5 x 10~ 7 s. As mentioned above, the decay probability 
is highest at I ~ D, thus the LLP with lifetime around 
t - 10- 6 (V20/fcTO)(m v /lGeV)(£:/10 2 Gey)- 1 s will in- 
duce more di-muon events in detector. 



IV. DISCUSSIONS AND CONCLUSIONS 

In this paper, we investigated the possibility of search- 
ing long-lived particle (LLP) produced by high energy 
cosmic ray colliding with atmosphere. The LLP may 
penetrate the atmosphere and decay into a pair of muons 
near/in the neutrino telescope. Such muons can be 
treated as the detectable signal. This study is moti- 
vated by recent cosmic electron/positron observations 
by PAMELA/ ATIC/Fermi. The new data suggests 
new source of electron/positron which may come from 
O(TeV) dark matter. In order to understand the dark 
matter thermal history, new light O(GeV) particles have 
been proposed. It is quite natural to conjecture that dark 
sector is complicated. There are more than one light par- 
ticles in dark sector, for example the dark gauge boson 
A' and associated dark Higgs boson hi . In this paper, 
we studied the scenario with A' heavier than h' and h! is 
treated as LLP. 

We have studied the LLP production processes and 
found that the promising process is single A' production 
qq — > A' , and A' subsequently decays into hi rather than 
into SM particles. Our numerical calculations show that 
for A' with mass of lGeV, the production rate can reach 
10 4 km~ 2 s~ 1 sr~ 1 and the final di-muon rate from h! is 
serval tens for some favorable parameter region. We have 
assumed the n ~ 10~ 3 , a' = a and the branching ratio 
of A' — > hi — > is 1. For different parameter se- 

lections, our results need to be multiplied by a factor of 
(k /10~ 3 ) 2 Br A>-+h' Br h , ^ ^+ ^- . Here it is worth remark- 
ing that a complete analysis for LLP production needs 
QCD correction to production process and simulations of 
cosmic ray shower included secondary hadrons and nu- 
cleons. The LLP could also be produced from secondary 
meson decay directly if allowed by kinematics [28| . These 
elements will increase LLP production rate efficiently. 

We simulated the signal di-muon events and calcu- 
lated the rate as functions of di-muon energy and the 
separation angle between two muons. Our numerical 
results showed that the large volume neutrino detector 
IceCube is suitable to detect LLP with lifetime about 
10~ 8 ~ 10~ 4 s. It is worth to mention that such parame- 
ter region could be compatible with the constraints^ from 
fixed-target experiments. For example, the Ref. [f| has 
reported the constraints on LLP decay length cr between 
1 cm and 10 8 cm by CHARM experiment result (29| . 

The scenario proposed in this paper could be cross- 
checked (tested) at the low energy e + e~ collider and/or 
large hadron collider. At the collider, the LLP will es- 
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cape from the detector and act as the missing energy. 
For example, the well promising process is the jA' asso- 
ciated production [3(| ■ The j + fl T signal can be isolated 
from SM irreducible background 7Z — > ^vv. If the LLP 
associated production with other light bosons which de- 
cay into charged leptons, the multi- e + e~ / ' fi + \i~ +f T is 
a cleaner signal. Moreover, the interaction between SM 
and dark sector may be induced via the mixing of Higgs 
fields. It implies that SM Higgs may decay into LLP. 
Such possible invisible decay modes can even change our 



search strategies of SM Higgs boson. 
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